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Worm-like filaments that are propelled homogeneously along their tangent vector are studied 
by Brownian dynamics simulations. Systems in two dimensions are investigated, corresponding to 
filaments adsorbed to interfaces or surfaces. A large parameter space covering weak and strong 
propulsion, as well as flexible and stiff filaments is explored. For strongly propelled and flexible 
filaments, the free-swimming filaments spontaneously form stable spirals. The propulsion force 
has a strong impact on dynamic properties, such as the rotational and translational mean square 
displacement and the rate of conformational sampling. In particular, when the active self-propulsion 
dominates thermal diffusion, but is too weak for spiral formation, the rotational diffusion coefficient 
has an activity-induced contribution given by Uc/^p, where Vc is the contour velocity and 5,p the 
persistence length. In contrast, structural properties are hardly affected by the activity of the 
system, as long as no spirals form. The model mimics common features of biological systems, such 
as microtubules and actin filaments on motility assays or slender bacteria, and artificially designed 
microswimmers. 


I. INTRODUCTION 

Its importance in biology and its enormous potential 
impact in technical applications makes active soft matter 
a field of rapidly growing interest and progress.Flexi¬ 
ble slender bodies are of particular importance. The ma¬ 
jority of natural swimmers propel themselves using flex¬ 
ible, hair-like structures like cilia and flagella.® Another 
important example are actin filaments and microtubules, 
major constituents of the cytoskeleton, whose capability 
to buckle decisively controls the mechanical properties 
of the cell body.® Flexibility is the crucial ingredient for 
the formation of small-scale spiral^ and possibly also 
for large-scale swirl^of microtubules on motility assays. 
Even the structure of slender bacteria can be dominated 
by their flexibility.® Elextrohydrodynamic convection can 
propel colloid chains because they are flexible,® just as 
the swimming mechanism of assembled magnetic beads 
in an oscillating external magnetic field is possible be¬ 
cause of the swimmer’s flexibility.® Elexibility is of course 
also the feature that allows for the instabilities leadi ng to 
cilia-like beating in artificially bundled microtubules.^^^^ 

Despite its importance, the number of theoretical stud¬ 
ies of active agents that incorporate flexibility is still 
relatively small, and can roughly be subdivided into 
works that focus on buckling phenomena and on free- 
swimming agents. Symmetry breaking instabilities lead¬ 
ing to rotation and beating motion of active filaments 
on motility assays can be described with a phenomeno¬ 
logical ordinary differential equation for the filaments.^ 
The propulsion force of motor proteins has been pre¬ 
dicted based on a Langevin model for buckled, ro- 
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tating actin filaments and microtubules .^1^ Numerical 
studies with Lattice-Boltzmann simulations and Brow¬ 
nian or multi-particle collision dynamics have demon¬ 
strated that clamped of pinned filaments composed of 
stresslets o r pro pelled beads can show cilia-like beating 
or rotation.EI^ 

The behaviour of free-swimming actin filaments on 
motility assays was reproduced in early numerical stud¬ 
ies using the Langevin equation.^ However, it was only 
recently that theoretical study of flexible, active fila¬ 
ments that can move freely has received significant atten¬ 
tion. Lattice-Boltzmann simulations reveal that sponta¬ 
neous symmetry breaking in chains of stresslets can lead 
to rotational or translational filament motion.l^ Brown¬ 
ian dynamics simulations of short self-propelled filaments 
suggest that different types of motion occur for single 
filamentJI^ and that spontaneous rotational motion can 
arise for pairs of filaments.^ A combination of Brow¬ 
nian dynamics simulations and analytic theory shows 
that shot noise in worm-like filaments leads to tempo¬ 
ral superdiffusive filament movement and faster-decaying 
tangent-tangent correlation functions.!^ Einally, chains 
of active colloids connected by springs have the same 
Elory exponent but a different prefactor of the scaling 
law compared to chains of passive colloids, as shown re¬ 
cently both analytically for beads without volume exclu¬ 
sion and numerically with Brownian dynamics simula¬ 
tions for beads with volume exclusion.^ 

The free-swimming behaviour of a worm-like filament 
that is tangentially propelled with a homogeneous force 
is still unexplored and is the subject of this work. The 
model is introduced in Section mi Results for the struc¬ 
tural and dynamic properties over a wider range or 
propulsion forces and filament flexibilities are presented 
in Section mg We find that the filament can sponta¬ 
neously form spirals, which is the mechanism that domi¬ 
nates the behaviour for large propulsion forces. The rele- 
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vance of our observations for natural and artificial active 
agents is discussed in Section [IV| We present our conclu¬ 
sions in Section \V\ 


II. MODEL AND METHODS 



We study a single, active, worm-like filament, which is 
modelled as a sequence of A +1 beads connected via stiff 
springs. The over damped equation of motion is given by 


FIG. 1. Filament model: Beads are connected via stiff 
springs. The active force acts tangentially along all bonds. 
Colour gradient indicates the force direction. 


7ri = -Vf/ + F«T+F«, (1) 

where are the coordinates of bead i, 7 is the friction 
coefficient, U is the configurational energy, is the 

thermal noise force, and Fp ^ is the active force that drives 
the system out of equilibrium. The configurational po¬ 
tential energy 


Without propulsion force, Fp ^ = 0, the model matches 
the well- known worm-like chain model for semi-flexible 
polymers) ^^ * ^^ * For active filaments, we use a force per 
unit length fp that acts tangentially along all bonds, i.e., 

N 


U — td)ond T b^angle T b^EV (2) 

is composed of a bond contribution between neighbouring 
beads 


b^bond 
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a bending energy 
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and an excluded volume term modelled with repulsive 
Lennard-Jones interactions 
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where Vij = — vj is the vector between the position 

of the beads i and j, ks is the spring constant for the 
bond potential, pq is the equilibrium bond length, hi is 
the bending rigidity, and e and a are the characteristic 
volume-exclusion energy and effective filament diameter 
(bead size). 

The drag force is the velocity of each bead times the 
friction coefficient 7 . The thermal force is modelled 

as white noise with zero mean and variance 2kBTj/At 
as described in Ref. [22l Note that hydrodynamic inter¬ 
actions (HI) are not included in our model. The model 
is thus in particular valid for (i) neutral swimmers, for 
which HI are known to be of minor importance,^^^J ^ (ii) 
swimmers near a wall, where HI is of less importancej ^^ * ^^ * 
and (hi) microorganisms t hat g lide on a surface, such as 
nematodes like C. elegansW^ 


as illustrated in Fig. The force along each bond is 
distributed equally onto both adjacent beads. 

We consider systems with parameters chosen such that 
(i) ks is sufficiently large that the bond length is approxi¬ 
mately constant ro, that (ii) the local filament curvature 
is low such that the bead discretization does not vio¬ 
late the worm-like polymer description, and that (iii) the 
thickness of the chain has negligible impact on the re¬ 
sults. When these requirements are met, the system is 
fully characterized by two dimensionless numbers. 


^p/L = 


K 

UbTL' 


Dt UbT ’ 


(8) 

(9) 


where L = Ntq and are the length and persistence 
length of the chain, respectively. ^p/L is a measure for 
the bending rigidity of the filament. The Peclet number 
Pe is the ratio of convective to diffusive transport and 
measures the degree of activity. For its definition, we use 
that the filament has a contour velocity Vc = fp/^h 
that the translational diffusion coefficient Dt = kpT/^iL^ 
where we have introduced the friction per unit length 

7i = + P/L. 

The ratio of these numbers 


^ = PeLlCp = ^, 


( 10 ) 


which we call the flexure number, provides a ratio of ac¬ 
tivity to bending rigidity. Previous studies showed that 
this num ber is decisive for buckling instabilities of active 
filamentit will be shown below that this is also a 
determining quantity for spiral stability and rotational 
diffusion. 

Simulations were performed in two dimensions, where 
volume exclusion interactions have major importance. 
Equations of motions were integrated using an Euler 
scheme. Simulation parameters and results are reported 
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FIG. 2. Trajectories of the center of mass of the filament (red) and filament configurations from selected snapshots (grayscale, 
leading tip black). Arrows point in the direction of movement. ^pjL — 0.2 (with N — 100) in all plots; Pe increases from left 
to right. At Pe — 200 (left), there is not sign of spiral formation. At Pe — 1000 (middle) spirals form occasionally, but the 
overall behaviour is dominated by an elongated chain. At Pe — 5000 (right), the spiral state is predominant. The chain has 
a directed motion in the elongated state. In the spiral state, the translational motion is almost purely diffusive. This leads to 
separated clusters in the trajectories for simulations in the strong spiral regime, visible for example in the upper left of the right 
image. The length of the depicted trajectories corresponds to approximately O.lSr (left), O.lSr (middle), and O.Gr (right). 


in dimensionless form, where length are measured in units 
of the filament length L, energies in units of the thermal 
energy /c^T, and time in units of the characteristic time 
for the filament to diffuse its own body length 

r = L^-iilAkeT. (11) 

In our simulations we used ks = 4000 fesT/rg, tq = a = 
L/N^ and e = ksT if not stated otherwise. A large pa¬ 
rameter space for Pe and ^p/L was explored by varying 
/p, N, and k. N was varied in the range from 25 to 200 
from the highest to the lowest ip/L. Almost all simu¬ 
lations were run for more than 5r. An initial period of 
the simulation output is discarded in the analysis. The 
timestep At was adjusted to the remaining settings to 
ensure stable simulations. Unless explicitly mentioned, 
results refer to simulations that were started with a per¬ 
fectly straight conformation. 

All simulations were performed using the LAMMPS 
molecular simulation packag^^ with in-house modifi¬ 
cations to describe the angle potential, the propulsion 
forces, and to solve the overdamped equations of motion. 

III. RESULTS 

The characteristic filament behaviour depends on its 
bending rigidity and activity and can be divided into 
three regimes (see Fig. [^. At low Pe or high ^pjL^ 
the “polymer regime”, the active filament structurally 
resembles the passive filament with Pe = 0. The main 
difference compared to the passive filament is that the 
active force drives the filament along its contour, lead¬ 
ing to a directed translational motion — we name this 
characteristic movement “railway motion”. At high Pe 
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(c) Spiral break-up by widening 


FIG. 3. Spiral formation and break-up mechanisms. Numbers 
in each panel provide the elapsed time. Leading tip is black. 
N — 100 in each case. (a,b) ^pjL — 0.4, Pe — 10 000, 5 = 
25 000. (c) ipjL = 1.0, Pe = 900, 5 = 900. 


and low ^pjL^ the filament spontaneously winds up to 
a spiral. The “spiral state” is characterized by ballistic 
rotation but only diffusive translation. Spirals can spon¬ 
taneously break up. Their lifetime determines whether 
spiral formation has a minor impact on the overall fila¬ 
ment behaviour, the “weak spiral regime” at intermediate 
Pe, or whether spirals are dominating, the “strong spi¬ 
ral regime” at large Pe. Because spiral formation has a 
major impact on both the structure and the dynamics, 
features related to spiral formation are addressed first. 
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FIG. 4. (a) Definition of the spiral number 5. The measure can effectively distinguish between elongated (red), weakly wound¬ 

up (blue), and strongly wound-up (orange) conformations, (b) Evolution of the spiral number 5 for the examples given in Fig.|^ 
with ^p/L = 0.2 (N = 100). Lines centered around 5 = 0 display results for a — LjN; lines centered around 5 = 8 (shifted 
upwards for clarity) show results for a = ILjN. (c) Probability distributions of the absolute value of the spiral number p(|5|) 
for the same examples with a — LjN. (d) Kurtosis ^2 for selected filament rigidities. Circles: a = L/N^ triangles: a = ^LjN. 
(e) Phase diagram. Background colour: kurtosis. Blue circles: polymer regime. Cyan triangles: weak spiral regime. Light and 
dark red squares: strong spiral regime. For the dark red squares, spirals did not break up during the simulations once formed. 
Black lines are a guide to the eye. Green area: threshold for spiral stability against break-up by widening. Spirals above this 
threshold will unfold by widening, spirals below will not. Purple: parameter space that can be obtained by actin filaments on 
a myosin carpet at T = 300 K using parameters for fp and n from Ref. [12] 


Structural and dynamic properties of the elongated and 
spiral state are presented afterwards. 


III.l. Spiral Formation 

The processes that lead to the formation and break-up 
of spirals are depicted in Fig. Spontaneous spiral for¬ 
mation (cf. Fig. Hf) results from the leading tip colliding 
with a subsequent part of the chain. Volume exclusion 
then forces the tip to bent. By further forward move¬ 
ment, the chain winds to a spiral. Two spiral break-up 
mechanisms occurred in our simulations. The first is the 
thermally activated mechanism in Fig. [^. The leading 
tip of the wound-up chain spontaneously changes direc¬ 
tion and the spiral deforms. This break-up mechanism 
requires strong local bending and is therefore predomi¬ 
nant for small ^p. The second mechanism is spiral break¬ 
up by widening and is depicted in Fig. The bending 
potential widens the spiral until the leading tip looses 
contact to the filament end. This break-up mechanism is 
predominant when ^p is too large for spontaneous spiral 


break-up. Because high stiffness is also unfavourable for 
spiral formation, spiral break-up by widening was almost 
exclusively observed in simulations that started with a 
spiral configuration. 

To understand spiral formation more quantitatively, 
we introduce the spiral number 

s = (0(L) - 0(O))/27r, (12) 

where (j){s) is the bond orientation at position s along 
the contour of the filament, as measure for the instan¬ 
taneous chain configuration. The definition is illustrated 
for three sample structures in Fig. It effectively mea¬ 
sures how often the filament wraps around itself. The 
time evolution of 5 is depicted in Fig.|^ for the same Pe 
and ^p/L as in Fig. At Pe = 200, 5 is always close 
to zero. At Pe = 1000, 5 behaves similarly, except that 
peaks with larger values for |5| occur occasionally, i.e., 
when spirals with a short lifetime form. At Pe = 5000, 
extended plateaus develop at large |5|. The spirals are 
wound up stronger and have a much longer lifetime. 

Probability distributions p(|s|) are depicted in Fig.|^. 
For the simulation without spirals (Pe = 200), the his- 
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togram resembles the right half of a Gaussian distribu¬ 
tion. For the simulation in the weak spiral regime, p{5) is 
similar for low |5|, but also has a small peak at |5| ~ 2 — 3. 
For strong spiral formation at Pe = 5000, p{\5\) has only 
a small peak at low |5|, which corresponds to the elon¬ 
gated state, and a large peak at large |5|, the predomi¬ 
nating spiral state. It turns out that the different regimes 
can be well distinguished by the kurtosis 



where (...) denotes the ensemble average and a 5 is the 
standard deviation of 5 . Results for the kurtosis are 
shown in Fig. W for selected ^pfL. /32 ~ 3 in the poly¬ 
mer regime, as expected for Gaussian distributions. In 
the weak spiral regime, the small peak at larger values 
increases the numerator in Eq. (13) and has only a weak 
impact on cr^, leading to an increase of ^ 2 - When the 
spiral state is dominating, grows drastically, resulting 
in a much smaller kurtosis ^ 2 - Note that to reduce sta¬ 
tistical uncertainties we symmetrized the 5-distribution 
in the computation of P 2 by counting each measured |5| 
as +5 and —5 and only used data from the spiral state 
for simulations that do not show spiral break-up. 

With the kurtosis as measure to characterize spiral for¬ 
mation, a phase diagram can be constructed as depicted 
in Fig. I§. Low filament rigidity and high propulsion 
Pe is beneficial for spiral formation. In particular, for a 
fixed propulsion strength per unit length, any chain will 
form spirals if it is sufficiently long, because increasing 
the chain length without modifying any other parameter 
corresponds to moving to the lower right in the phase 
diagram. 

The dimensionless numbers ^p/L and Pe completely 
characterize the system if the filament diameter — or 
the filament aspect ratio — is of minor importance. This 
is the true in the entire polymer regime, where volume- 
exclusions interactions hardly come into play because of 
the elongated chain structure. For the spiral regimes, the 
aspect ratio has an impact on the structure of the spiral 
and does in this way influence the results. Which features 
of the spiral regime can be approximated well by the 
dimensionless numbers can be understood from the spiral 
formation and break-up mechanisms. The aspect ratio is 
hardly relevant for spiral formation and spiral break-up 
by widening, where the decisive moments are when the 
filament tip collides with subsequent parts of the chain, 
or when it looses contact to the chain end, respectively. 
That break-up by widening is characterized well by the 
dimensionless numbers is also confirmed by a series of 
simulations that we start from a spiral configuration in 
which we vary fp, /^, and kpT. It turns out that 
spirals will break up by widening if 


d < 1000 - 1500. 


(14) 


In contrast, spontaneous spiral break-up by a change 
of orientation of the leading tip is dependent on a strong 


local curvature close to the tip and the structure of the 
spiral, which in turn is dependent on the filament diam¬ 
eter. The dimensionless description does therefore not 
provide a full characterization of the strong spiral regime, 
where spontaneous spiral break-up is the only mechanism 
to escape the spiral state. This is also confirmed by re¬ 
sults for spirals that never broke up (cf. dark red squares 
in Fig. results for ^p/L = 0.2 and ^p/L = 0.14 
show non-mo notonic behaviour in the direction of ^p/L. 
This is a result of a combination of that the dimension¬ 
less description is only partially valid in this regime and 
that we chose N = 200 for ^p/L < 0.2 but N = 100 
for 0.2 > ^p/L > 2.0 in our simulations, i.e., the aspect 
ratio L/a is halved in our simulations for filaments with 
Cp/L < 0.2. 

Finally, the bead discretization with the chosen pa¬ 
rameters favours a staggered ar range ment of beads of 
contacting parts of the filamentwhich implies an 
effective sliding friction between these parts. To study 
the importance of this effect, we increase the diameter 
of the beads at fixed bond length so that neighboring 
beads are heavily overlapping, which leads to a strongly 
smoothened interaction potential. Results for an in¬ 
creased diameter a = 2LfN are shown in Fig. & and d. 
We find that the spiral formation frequency is hardly af¬ 
fected by smoothening the filament surface. In contrast, 
spontaneous spiral break-up is largely alleviated for the 
smoother filament, leading to decreased spiral life-time, 
as can be seen from the evolution of 5 for Pe = 5000 
in Fig. Smoother filaments thus show a qualita¬ 

tively similar phase behaviour with slightly moved phase 
boundaries. 


III.2. Structural Properties 


The structural properties of the filament conformations 
can be best understood from the end-to-end vector rg, as 
depicted in Fig. As long as no spirals form, simula¬ 
tion results are in good agreement with the Kratky-Porod 
model (valid for worm-like, non-active poly mers without 
volume-exclusion interactions) that predict^^^^ 


in two dimensions. At low ^p/T, volume-exclusion in¬ 
teractions lead to slight deviations between the Kratky- 
Porod model and the simulation results. Strong devia¬ 
tions between the Kratky-Porod model and simulation 
results only occur in the strong spiral region in the phase 
diagram. The same trend was observed for the tangent- 
tangent correlation function, the radius of gyration, and 
the static structure factor, but is not reported here to 
avoid unnecessary repetition. 
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FIG. 5. Symbols: Mean end-to-end distance \/(r^ over Pe 
for different values of ^pjL. Solid lines: y/ (r^) as predicted 
from the Kratky-Porod model. The symbol shape indicates 
the region in the phase diagram: circle: polymer regime; tri¬ 
angle: weak spiral regime, squares: strong spiral regime. N 
varies from 25 to 200 from large to small ^pjL. 


III.3. Dynamic Properties 


The characteristic filament motion can be understood 
from the mean square displacement (MSD) of a bead jpi 
relative to bead as shown in Fig. For comparison, 
the MSD of the reference bead j of a passive filament 
is also shown. Note that the displacement of this bead 
is subdiffusive at the short lag times shown here.l^ Dis¬ 
placement functions of the propelled beads are always 
larger than in the passive case. The curves show three 
distinct regimes. At small lag times, the MSDs of active 
filaments display plateaus due to the average distance 
of the two beads j + i and j along the filament. At 
large lag times, the increased motion caused by activity 
effects the MSDs of the propelled beads to grow more 
rapidly than that of the passive bead. The relevant part 
of the MSD that shows that the characteristic filament 
motion is movement along its contour is at intermediate 
lag times, where the MSDs of the propelled beads pass 
through minima that touch the reference MSD for ther¬ 
mal motion. At that lag time, the bead j i has moved 
approximately to the position of bead j at zero lag time. 
The beads have moved along the chain contour, similar to 
the movement of a train on a railway. The deviation from 
the exact starting position of bead j exactly matches the 
thermal motion, which results in the MSDs of the pro¬ 
pelled beads touching the MSD of the passive filament. 
Thus, the characteristic movement of the filament is mo¬ 
tion along its contour superimposed with thermal noise, 
as depicted in Fig. IT- Note that ^p/L = 0.3 was selected 
in Fig. [^because this corresponds to a rather flexible fil¬ 
ament, where stronger deviations from the characteristic 
railway motion might be expected. This type of motion 
was observed in all simulations in the polymer regime 
and weak spiral regime. In the strong spiral regime, the 
MSDs of the propelled beads even fall below the reference 
line for purely thermal motion. 



FIG. 6. Mean square displacement of bead 25+ z with respect 
to bead j = 25 {N = 100). ^pjL — 0.3 for all lines. Black 
line: Pe — 0, coloured lines: Pe — 1000 (weak spiral regime). 
Black line describes purely diffusive motion of the leading 
bead j. 




FIG. 7. (a) Series of snapshots of a hlament with jL — 0.3 
(with N — 100) and Pe — 1000. The chain moves along 
its contour superimposed with thermally activated motion. 
Golouring is to improve the distinctness of the chains, (b) Ide¬ 
alized railway motion in the absence of diffusion. The hlament 
(thick black line) moves with velocity Vc along the contour of 
an inhnite chain with same ^p (gray line). so(t) runs along 
the contour of the inhnite chain and marks the end point of 
the hlament. 


The rotational diffusion can be accessed from the ori¬ 
entation of the end-to-end vector 0. Its mean square 
rotation (MSR) is given in Fig. |^. Note that complete 
rotations around the axis are accounted for in our com¬ 
putations. 0 {t) can therefore be much larger than 27r. 
In both the spiral and the elongated state, there is a 
regime at short lag times in which the MSR is domi¬ 
nated by the internal hlament hexibility. For the spiral 
state, this regime is followed by a ballistic regime with 
MSRoc For simulations in which the spirals break up 
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(b) 


FIG. 8. (a) Mean square rotation of the end-to-end vector; 
the red line is from the strong spiral regime, the blue from 
the polymer regime {N = 100). (b) Rotational diffusion co¬ 
efficient Dr as a function of the flexure number 5. Symbol 
shape indicates the region in the phase diagram: circle: poly¬ 
mer regime; triangle: weak spiral regime. Gray fines are pre¬ 
diction with Eq. (24) for different Dr,p; for the lowermost fine 
Drp = 9/4r“^. N varies from 25 to 200 from large to small 
^p/L. 


spontaneously, a subsequent regime at high times with 
MSRoc t is expected but could not be detected in our 
simulations because of finite simulation time and strong 
noise at large lag times in the MSR. For the elongated 
state, the regime dominated by internal flexibility is fol¬ 
lowed by a diffusive regime with MSRoc t. 

The rotational diffusion coefficient Dp can be extracted 
by fitting MSR= 2 Drt to the regime of the MSR with 
gradient unity on a double-log scale. Measured Dp are 
given in Fig. as a function of the flexure number 
The diffusion coefficients collapse to a single curve, which 
has a plateau at low ^ and then grows linearly. Strong 
deviations from this trend are only observed at high ^ 
when the filament is in the weak spiral regime, and at 
^ = Pe = 0 for flexible filaments, where strong deviations 
from a rod-like shape increase Dp. 

The rotational diffusion coefficient Dp can be predicted 
from the characteristic railway motion in Fig. [^and the 
relation of the rotational diffusion coefficient to the au¬ 


tocorrelation function of the end-to-end tangent vector 

te 


{te(i)-te(0))=e-^''*, (16) 


which is valid for lag times t that are sufficiently large 
such that variation of te is not dominated by non- 
diffusive behaviour at early lag times caused by the fila¬ 
ment flexibility (cf. Fig. [^). With 


where t(s,t) is the tangent vector at position s of the 
filament at time t, the left hand side of Eq. (16) becomes 


{te{t) ■ te{0)) = ^ J ds' j ds{t{s,t) ■t{sl,0)) , 

(18) 

where the order of summations has been changed to ar¬ 
rive at the right-hand side of Eq. Q. As a representa¬ 
tion of the characteristic railway motion (cf. Fig. [^), we 
write 


t(5, t) = t{s + Vet, 0). (19) 


Note that this equation disregards the passive equilib¬ 
rium rotation Dp^p. With Eq. (19) and the expression 


for the ta ngent -tangent correlation function of worm-like 
polymers,!^®^ the integrand in Eq. (18) becomes 

(t(s, t) • t{s', 0)) = (t(s + Vet, 0) • t{s', 0)) 

= exp[-(s + Vet - s')/^p]. (20) 


Integrating Eq. (18) provides 


(te(l) • te(0)) = (2 - - e--^/^^)) . 

( 21 ) 

A second order Taylor expansion in (small) L/^p then 
gives 


(te(t) • te(0)) = ex.p[-Vct/^p], (22) 

so that a comparison with Eq. finally yields the 

activity-induced rotational diffusion 

Dp^a='^c/^P- (23) 

Note that Vc/^p = 5^/4r. Assuming that uncorrelated 

activity-induced and thermal rotation Dp^a and Dp^p con¬ 
tribute to the overall rotation, we write 

Dp = Dp^p + Dp^ai (24) 

where Dpp depends on ^p/L and has the lower bound 
Dp^p = (9/4) for rod-like filament s.l^ As can be seen 
from Fig. [^, Eq. matches the simulated rotational 
diffusion coefficient accurately. 

The characteristics of the center-of-mass MSD are 
shown in Fig. For the polymer regime, the typical 



























FIG. 9. Mean square displacement of the center of mass. 
Dashed coloured lines : si mulation results. Gray solid lines: 
predictions using Eq. ( [^ . (a) Results for = 4.0 (polymer 
regime, N = 50) and different values of Pe increasing from 0 
to 100 000. vq is predicted from Eq. ( [27t , Dr from Eq. ( [24| 
with Dr,p — . (b) Results for ^pjL — 0.3 (spirals at 

large Pe, N — 100). vq from Eq. (27), Dr is determined from 
a ht to the measured MSD, because Eq. ( [^ is only valid in 
the polymer regime. Spiral formation leads to a decreased 
MSD. 


S-shape of subsequent short-time diffusive, intermediate¬ 
time ballistic, and long-time effective diffusive behaviour 
developspl stronger propulsion increases the MSD. An 
important difference compared to rigid bodies is that 
the transition time = 1 /Dr to long-time diffusive be¬ 
haviour is dependent on the propulsion strength. 

When spiral formation becomes important, the general 
trend of the MSD changes, as shown in Fig. for a 
flexible filament. In the polymer regime or weak spiral 
regime, increasing Pe leads to a larger displacement. In 
the strong spiral regime, however, the MSD decreases. 
For very stable spirals, the MSD is only weakly affected 
by the propulsion and almost matches the case of purely 
diffusive motion. 

The MSD for active point particles, spheres, or stiff 
rods is given 

{{rdt) - rc(0))^) = 4Dtt + 

{2vo/D'^)[Drt + exp{-Drt) - ^^5) 


FIG. 10. {si)S{q,t)/S{q,0) for ^p/L = 1.0, N = 100, q « 
Stt/L and different Pe. (b) rs{q) for q ^ Stt/L . Gircles 
correspond to the polymer regime, triangles to the weak spiral 
regime, and squares to the strong spiral regime in both plots. 


where Dt is the translational diffusion coefficient and uq 
is a ballistic velocity. It turns out that Eq. (25) can be 
used to describe the MSD for active filaments, when the 
three coefficients Dt^ uq, and Dr are chosen properly. 
The translational diffusion coefficient is Dt = L‘^/At = 
ksT/^iL. We predict the rotational diffusion coefficient 
Dr with Eq. ( [2^ . Finally the effective velocity can be 
expressed via 



(26) 


as a balance of the net external force \Fp\ with the total 
friction force ^iLvq. |Fp| can conveniently be expressed 
as the propulsive force per bond fp times the end-to-end 
vector, thus leading to 


no = 


Sy-JW) 


(27) 


As shown in Fig. using these correlations for the coef¬ 
ficients provides an accurate prediction of the MSD. 

The last item we address is the effect of propulsion on 


conformational sampling. Figure 10 r shows results for 
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the dynamic structure factor 

/ AT+l AT+l \ 

= ( tvTI ^ ^ exp{iq ■ [rj(0) - rj(i)]} \ 

\ ^ i=l j = l / 

(28) 

averaged over different directions of q. In the phase 
without spirals, 5'(g, t)/5'(g, 0) decays more rapidly with 
increasing Pe, indicating a faster change of conforma¬ 
tions with increasing propulsion. When spirals form, 
S{q,t)/S{q,0) is indepenent of Pe and larger than 
S{q, t)/S{q, 0) at Pe = 0, indicating a slow change of con¬ 
formations, which agrees with the observation of hardly 
any internal motion of the chain in this regime in our sim¬ 
ulation output. Note that for the strong spiral regime, 
the data is from simulations where spirals formed spon¬ 
taneously and did not break up. The depicted data is a 
result of averaging over the spiral states only. 

To better quantify the behaviour of S{q^ t), we compute 
the characteristic decay time of the dynamic structure 
factor 


Ts{q) 


tS{q, t)dt 
f^S{q,t)dt ■ 


(29) 


Results for rs{q) at g ^ Stt/P, a g-vector large enough to 
capture the behaviour of mainly the internal degrees of 
freedom, are given in Figjl^. rs{q) decays slowly at low 
Pe. At high Pe, ts decays inversely proportinal to Pe 
when no spirals form This is consistent with the picture 
that instantaneous conformations are essentially identical 
to those of passive filaments, but they are traversed with 
velocity Vc, corresponding to r (x Pe“^. In the strong 
spiral regime, ts is large and independent of Pe, which 
is a sign for that conformational changes are irrelevant 
and that rs is determined by the quasi-diffusive center of 
mass movement. Note that the measured ts at different 
^p/L collapse to a single line for both the polymer and 
the strong spiral regime. 


IV. DISCUSSION 

The spontaneous formation of spirals is the feature 
dominating the overall behaviour of self-propelled fila¬ 
ments, both for dynamic and structural properties. For¬ 
mation of spirals was previously observed for long, slen¬ 
der bacteria surrounded by short bacteria.l^ It was con¬ 
cluded that interaction with other active particles is a 
prerequisite for spiral formation. In contrast, the study 
at hand shows that spirals can form even for isolated fil¬ 
aments, as long as (i) the filament is sufficiently flexible, 
(ii) the propulsion is sufficiently strong, and (hi) excluded 
volume interactions force the tip of the filament to wind 
up. 

The first two conditions will be met automatically for 
any real system by choosing L sufficiently large and leav¬ 
ing all other parameters constant (leading to increased 
Pe and decreased ^pjL^ i.e., favouring spirals). Meeting 


the third condition can in general not be achieved so eas¬ 
ily. A free-swimming filament in three dimensions or a 
filament in two dimensions with low resistance of cross¬ 
ing its own body will not form spirals. This is also one 
reason why spiral formation has not yet been observed 
in more experimental studies. Agents that are similar to 
our model are actin filaments or microtubules on a motil¬ 
ity assay. The former have a high crossing probability,!^ 
formation of spirals is therefore not expected. The area 
enclosing the act in-filament parameter space in Fig. 
must thus be understood as that the regime where the 
flexibilities and propulsion strengths permit spiral for¬ 
mation can in principle be reached in real systems, and 
not so much as that sufficiently long actin filaments will 
form spirals. Microtubules on dyn ein c arpets, which have 
a much lower crossing probabihty,^^^ will possibly form 
spirals if they are grown to sufficient size. 

Overall, except for slender bacteria,!^ we are unaware 
of a microscopic example in which spiral formation was 
observed. Yet, the formation of spirals is a feature that 
deserves more attention in the future. First, formation 
of spirals is an extremely simple non-equilibrium phe¬ 
nomenon that, in contrast to many other phenomena of 
active matter, arises for a single self-propelled particle 
and cannot easily be mapped qualitatively to passive sys¬ 
tems in which activity is replaced by attractive forces. It 
can thus be used as a model phenomenon for the study 
of non-equilibrium thermodynamics. Second, our model 
is very simple; a realization in experiment seems possible 
within the near future. Finally, the formation of spi¬ 
rals leads to a sudden, strong change in structural and 
dynamic properties. The effect can thus potentially be 
used as a switch on the microscopic scale. 


V. SUMMARY AND OUTLOOK 

We report an extensive study for the behaviour of di¬ 
lute, self-propelled, worm-like filaments in two dimen¬ 
sions. The spontaneous formation and break-up of spi¬ 
rals is the feature that dominates the filament behaviour. 
Spiral formation is favoured by strong propulsion and 
low bending rigidity. Propulsion has a noticeable impact 
on structural properties only when spirals are dominat¬ 
ing. The Kratky-Porod modeP^ is therefore valid for 
filaments that are weakly propelled or have high bend¬ 
ing rigidity. When spiral formation becomes significant, 
structural properties change drastically. 

The characteristic filament motion is what we call the 
railway behaviour. The chain moves along its own con¬ 
tour superimposed with noise. With the understanding 
of the structural properties and the characteristic mo¬ 
tion, rotational diffusion and the center-of-mass mean 
square displacement can be predicted to high accuracy 
when no spirals form. In contrast to rigid bodies, propul¬ 
sion has an impact on the rotational diffusion coefficient. 
Finally, propulsion enhances conformational sampling in 
the regime without spirals. 
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An obvious next step is understanding the collective 
motion of such active filaments. We expect that our sin¬ 
gle filament results will help to understand the collective 
behaviour, which is nonetheless strongly influenced by 
the additional interactions. In particular collision with 
other constituents might enhance spiral formation and 


lead to swirl-like patterns. 
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